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ABSTRACT 

14 Although propane gas (CsHg) was first detected in the stratosphere of Titan by the 

15 Voyager IRIS infrared spectrometer in 1980, obtaining an accurate measurement of its 

16 abundance has proved difficult. All existing measurements have been made by modeling 

17 the 1^26 band at 748 cm~^: however different analyzes over time have yielded quite different 

18 results, and it also suffers from confusion with the strong nearby band of acetylene. 

19 In this paper we select large spectral averages of data from the Cassini Composite In- 

20 frared Spectrometer (CIRS) obtained in hmb- viewing mode at low latitudes (30°S-30°N), 

21 greatly increasing the path length and hence signal-to-noise ratio for optically thin trace 

22 species such as propane. By modeling and subtracting the emissions of other gas species, 

23 we demonstrate that at least six infrared bands of propane are detected by CIRS, in- 

24 eluding two not previously identified in Titan spectra. Using a new line list for the range 

25 1300-1400 cm~^, along with an existing GEISA list, we retrieve propane abundances from 

26 two bands at 748 and 1376 cm~^ At 748 cm~^ we retrieve 4.2 ±0.5 x 10~^ (l-cr error) at 2 

27 mbar, in good agreement with previous studies, although lack of hotbands in the present 

28 spectral atlas remains a problem. We also determine 5.7 ± 0.8 x 10^^ at 2 mbar from the 

29 1376 cm~^ band - a value that is probably affected by systematic errors including contin- 

30 uum gradients due to haze and also an imperfect model of the vq band of ethane. This 

31 study clearly shows for the first time the ubiquity of propane's emission bands across the 

32 thermal infrared spectrum of Titan, and points to an urgent need for further laboratory 

33 spectroscopy work, both to provide the line positions and intensities needed to model 

34 these bands, and also to further characterize haze spectral opacity. The present lack of 

35 accurate modeling capability for propane is an impediment not only for the measurement 

36 of propane itself, but also for the search for the emissions of new molecules in many spec- 

37 tral regions. 
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39 1 Introduction 



40 Propane (CsHg) was first detected in the atmosphere of Ti tan by the Voyager 1 IRIS 



198ll ). when excess emission at 



spectrometer during the 1980 encounter (iMaguire et al. . 
748, 922, 1053 and 1158 cm^^ was used to make the identification. Propane is an end- 
product of Titan's active photochemistry. This process begins when methane is dissociated 
in the upper atmosphere by UV photons and energetic particles to form hydrocarbon 
radicals and ions, which undergo a complex chain of reactions to form larger molecules 
and haze particles. The main pathway leading to C.-^Hs is the three-body ass ociation: CH3 



+ C2H5 + M ^ CsHg + M flYung et al. 



1984 



Wilson and Atreya 



20041 ). and most is 



48 ultimately removed by condensation in the lower stratosphere. Spectrally, propane is a 

49 very important molecule, as it has 27 predicted infrared vibrational modes of which at 

50 least 23 are active (see §4.1). 

51 Early attempts to measure the abundance of propane from the Voyager dataset suffered 

52 from a lack of reliable laboratory spectroscopic data, resulting in initial over-estimates 

53 follow ed by a downward trend of revised abundances until the mid-1990s. 



Maguire et al. 



54 (1198 ll ) used laboratory spectra to compute the band intensity of the 748 cm ^ band 

55 thereby inferring a uniformly-mixed abundance for Titan's stratospher e of 2 x 10~^, with 



56 an un certainty factor of three. The value was halved the next year by 



Kim and Caldwell 



19821). who modeled the same band to obtain a propane abundance of 1.2 x 10 ^. 



Coustenis et a. 



59 sorptions from 



1989a) a: 



Giver et al. 



s o used a band-model approach, based on newly measured ab- 



(119841 ) , and a full radiative transfer treatment to greatly revise 

60 the estimate for uniformly-mixed stratospheric propane downwards to (7 ± 4) x 10^^ at 

61 the equator. 

In 1991 ~9000 lin es representing the 748 cm~^ band were added to the GEISA atlas 



63 ( Husson et al. 



I992I ). with the data taken from unpublished measurements at NASA 
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GSFC by G. Bjoraker circa 1986. These allowed ICoustenis and BezardI (119951 ) to make a 
full line-by-line model of the blended C2H2 and CaHg emissions. For the first time, the 
latitude variation of propane was reported: a modest increase from 5.0 ± 1.7 x 10^'' at 
the equator to 1.2 ± 0.5 x 10^^ at 70°N (assuming uniform stratospheric abundances). 
Simultaneously, enhancements in many other minor gases were detected in the northern 
(winter) polar stratosphere: especially in the heavier, acetylene- derived gases (C4H2, C3H4, 
HC3N, C2N2), but also HON and C2H4. These enhancements have been modeled and 
explained dynamically, as due to downward advection of short-lived species at the pole(s) 



in a global circulation cell(s), as these m olecules tend to have 



altitudes in the absence of dynamics (ITeanby et al 
references therein). 



2009 



ligher abundanc e s at h igher 



Lebonnois et al 



2009 



and 



75 The launch of the Infrared Space Observatory (ISO) in 1995 provided a new tool 

76 for studying planetary atmospheres. In 1997, the Short Wavelength Spectrometer (SWS) 

77 was used to acquire spectra of Titan with a mean resolution of 0.4 cm~^ - an order of 

78 magnitude improvement over the 4.3 cm~^ of Voyag er IRIS, although spatia lly unresolved. 



Coustenis et al. 



(l2003l ). and included 



79 Low-latitude gaseous abundances were reported by 

80 a measurement of gcsHs = 2 ± 1 x 10"'' via the same (748 cm"^) band as Voyager IRIS. 

81 The reason for this slightly lower value is unclear, although real seasonal change might 

82 be a factor, as well as the higher spectral resolution. 

83 Even this improved spectral resolution was insufficient to actually separate the emis- 

84 sion lines of propane at 748 cm~^ from those of the R-branch lines of the 729 cm~^ (z/5) 

85 band of C2H2, and weaker contributions from HON and C';iH R. The first study to resolve 

86 the CaHg emission features from the other gases was made by 

87 the R~ 10^ TEXES Echelle Spectrograph attached to the NASA Infrared Telescope Facil- 

88 ity (IRTF) to make high spectral resolution, disk averaged observations of Titan, resulting 

89 in a stratospheric abundance estimate of (6.2±1.2) x 10"'' at 90-250 km (13.0-0.24 mbar). 



Roe et al 



fl2003h . who used 
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90 in agreement with the low-latitude IRIS values. It is important to note that these authors 

91 compared four line lists for the 748 cm~^ propane band then available, noting significant 

92 discrepancies between all the lists and actual propane spectral data. Only the absorption 

93 coefficients that they derived empirically from laboratory data provided a good fit to the 

94 planetary data. In addition, this was the first study of Titan's propane to use the z/26 

95 designation for this band, rather than the z/21 designation that was commonly used by 

96 earlier authors. We discuss the dichotomy of propane band designations further in section 

97 mi 

98 In 2004 the Cassini/Huygens spacecraft entered orbit around Saturn, opening a new 

99 era in Titan research. The Composite Infrared Spectrometer (CIRS) instrument onboard 

100 Cassini is a successor to the Voyager IRIS, although with many improvements including 

101 to spectral range and resolution, and spatial resolution (see Section [2TT]) . Furthermore, 

102 the great variety of Titan flyby encounter ranges and relative orbital inclinations, com- 

103 bined with the much greater available observing time versus the Voyager encounters, 

104 permits Cassini CIRS to make a very diverse range of scientific investiga tions. These in 



105 elude low-spectral resolution North-So uth surface temperature mapping (jjennings et al. 



20091 ) and limb temperature mapping (lAchterberg et al 



20081 ). through to high spectra. 



107 resolution composition mea s urements on both th e disk and the limb (ICoustenis et al. 



2007 



Vinatier et al. 



2007a 



Teanby et al. 



2008bl . and references therein), with vertical 



109 resolution of better than one scale height (40-50 km). 

110 So far, two analyses of CIRS Titan spectra have considered propane. 



Coustenis et al. 



111 (120071 ) modeled zonally-averaged sets of mid-infrared spectra (600-1400 cm~^) from 70°S- 

112 70°N to measure the latitude variations in hydrocarbons, nitriles and CO2, including the 

113 748 cm~^ band of propane for which spectral line data is available in GEISA. CsHg was 

114 found to exhibi t a very slight increas e from the south to the north, from ~ 5 x 10^'' 



to 



8 X 10- 



Vinatier et al.l (l2007al ) modeled two high spectral resolution (0.5 cm ^] 
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116 observation sets of Titan's limb, from the Tb and T3 flybys and directed at 15°S and 

117 80°N respectively. CsHg was found to display an increasing abundance with altitude in 

118 the range ~150-350 km (~ 2.0-0.01 mbar) at both latitudes, in accordance with the 

119 conventional understanding of photochemical formation in the upper atmosphere and 

120 downward advecti qn to the lower stratospher e where removal through condensation occurs 



121 (89% according to 



Wilson and Atreya 



20M). 



122 In this paper we have two objectives. First, we show clearly the very prolific extent 

123 of the multiple propane bands that can be detected by CIRS across the mid-infrared. 

124 At least six bands are visible in limb spectra after the modeling and removal of other 

125 stronger gas emissions. Second, we use a new laboratory line list for propane in the 

126 range 1300-1400 cm'^ Jpiaud et all . boOll ). along with the GEISA list for the 748 cm'^ 

127 band to make independent measurements of the propane abundance, and show that the 

128 results are probably compatible when systematic errors are accounted for. Finally, we 

129 draw conclusions about the current status of propane retrievals in the thermal infrared. 



130 2 Data Acquisition 

131 In this paper we make use of two spectral datasets. The first is a set of Titan limb 

132 spectral data, acquired by Cassini CIRS between 2004 and 2008. The second is a set 

133 of laboratory absorption spectra of room temperature propane gas acquired at Pacific 

134 Northwest National Laboratory in 2002. Both datasets are now described. 

135 2. 1 Cassini CIRS Limb Observations 



136 The Composite Infrared Spectrome ter (CIRS) carried on-board Cassini is a dual 

137 Fourier Transform interferometer design (IKunde et al.l . 



1996 



Flasar et al. 



200J). A shared 
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138 telescope and foreoptics feed incoming radiation simultaneously into two interferometers: 

139 a Martin-Puplett (polarizing) type optimized for far-infrared rays (10-600 cm~^), and 

140 a conventional Michelson design for the mid-infrared (600-1400 cm^^). The recombined 

141 beam of the far-infrared interferometer is focused onto a single circular bolometer detec- 

142 tor, known as focal plane 1 (FPl), while mid-infrared detection capability is provided by 

143 two 1 X 10 HgCdTe arrays: focal plane 3 (FP3, 600-1100 cm"^) and focal plane 4 (FP4, 

144 nominal range 1100-1400 cm~^, see also Appendix A). Both interferometers also share a 

145 common scan mirror mechanism and reference laser. The spectral resolution is determined 

146 by the commanded scan time, and is variable from 15.5 to 0.5 cm~^. 

147 At 5-9 hours from closest approach to Titan (1.0-1.8x10^ km range), the CIRS in- 

148 strument makes observations of the atmospheric limb, placing the mid-infrared arrays 

149 perpendicular to the disk edge, and obtaining spectral samples at ~30-50 km resolution, 

150 equivalent to approximately one atmospheric scale height. These observations come in two 

151 types: (i) temperature maps (shorter dwell, 15.5 cm~^ spectral resolution, repositioning 

152 at multiple latitudes) and (ii) composition integrations (longer dwell, 0.5 cm~^ resolution, 

153 single latitude). 

154 For the purpose of identifying the weak emission bands, we have created large averages 

155 of low-latitude spectra for FP3 and FP4 from the high spectral resolution limb integrations 

156 (type (ii) above). Note that we do not attempt to explore the vertical and latitudinal 



157 variation o: 



the propane abundaiice from the CIRS data s et, w. 



158 elsewhere (ICoustenis et al 



20071 : 



Vinatier et al. 



lich has been documented 



2007a 



20091 ) . but rather to maximise 



159 signal-to-noise to emphasize and detect very weak emission features while still maintaining 

160 reasonable homogeneity of the dataset. 

161 Therefore we chose to average across all flybys from the prime mission (June 2004 to 

162 June 2008) and latitudes in the range 30°S-30°N, in the knowledge that Titan's strato- 
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163 sphere is slowly 



Achterberg et al. 



varying; in this temporal a nd spatial interval (jCoustenis et all 12007 : 



2008 



Teanby et al. 



2008al ). However, Titan's atmosphere varies much 



165 more rapidly in the vertical direction, and therefore we restricted our averages to 100- 

166 150 km (pixel centers) where the signal-to-noise ratio was found to be the most favorable. 

167 The spacecraft distance range was also constrained to be less than 2 x 10^ km. See Table 

168 I for additional details. 



[TABLE I] 



170 2.2 PNNL Laboratory Spectroscopy 

171 For comparison with the CIRS spectra, Dr Steven Sharpe at Pacific Northwest National 

172 Laboratories (PNNL), has shared with us room temperature (296 K) absorption spectra 

173 of propane gas, taken as part of a large-scale survey of molecular infrared spectra. These 

174 consist of 6 individual co-added spectra of 98% pure propane gas measured with a Bruker- 

175 66V FTIR at 0.112 cm~^ unapodized resolution, covering the FP3 and FP4 part of the 



176 CIRS spectral 



177 the literature (ISharpe et al 



range (600- 



400 cm ^). Full experimental details have been published in 



2004h 



178 We have convolved the PNNL spectra with a Hamming kernel h{Ai)) (instrumental 

179 line shape) of FWHM 0.48 cm~^ for approximate comparison with the CIRS spectra, 

180 which originate at much lower temperatures (~ 160 K) in Titan's stratosphere: 



a 1.08 - OMa^ [Ai>l^ sine (27raAi>) 

■■■ "f^^' - T^l^' « 

182 where h is the weighting factor at each position Az/ = (z/ — z/q) relative to a given 

183 wavenumber uq in the original spectrum, and a is an effective path difference: a = 

184 1.815/(2 X 0.48) = 1.89 cm that compensates the Nyquist path difference at the de- 
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185 sired resolution (1/(2 x 0.48) cm) for the FWHM of the Hamming function in Nyquist 

186 units (1.815). Note that h is normahzed so that: / /i(Az/) d{Ai)) = 1.0 over the desired 

187 intervah we used ±10 samples relative to each wavenumber (21-pt kernel). 



188 3 Data Analysis 

189 3.1 Model Atmosphere 



Our model atmosphere consists of 100 layers equally spaced in \og{p) from the sur- 
face (1.45 bar) to 700 km (4 x 10"^ bar.) The initial temperature-pres sure profile is a 



20051 ) (covering from 



192 smoothed version of the Huygens HASI results ( iFulchignoni et al. . 

193 the surface to 7 x 10~^^ bar at 10°S, 168°E) with heights calculated from hydrostatic 

194 equilibrium. All gases known to be relevant for the 700-1400 cm"^ range of Titan's in- 

195 frared spectrum are included, namely: N2, CH4, H2, HON, C2H2, C2H6, C2H4 and CaHg. 

196 In addition, the isotopologues ^'^C^^CH2, ^^CH4, CH3D and ^^CHsD are included as sepa- 



197 rate spe cies. The initia 



198 results (Niemann et al 



abundances and isotopic ratios mostly follow the Huygens GCMS 



199 itudes fICoustenis et al 



2005h and previous measurements by CIRS pertaining to low lat- 
20071 ). In particular, H2 and C2H4 are uniformly mixed, while 
CH4 follows the GCMS profile, N2 is defined as = 1.0 — qcHA, and the other gas 
primary isotopologues, with the exception of C2H2, have constant abundances in the 
middle atmosphere, and follow condensation profiles towards the tropopause. The C2H2 
VMR, which is know n to have a gradient at low-latitudes, follows the 10°S profile of 



Vinatier et al 



(l2007al ). Finally, the minor isotopologue species have initial profiles scaled 



205 from the primary isotopologue by molecu lar isotopic ratios previously determined by CIRS 



(Bezard et al.. 


2007; 


Nixon et al.. 


2008a) 



2008al ). Fig. [T] shows the initial gas profiles. 



Fig. 1. Appears Here. 
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207 Finally, our model atmosphere also includes a single aerosol component with c o nstan t 

208 mass mixing ratio of particles above the tropopause (45 km) as described in iNixoru (119981 ). 



209 3.2 Forward Spectral Model 



The forward model (NEMESIS, llrwin et al.l . |2008| ) compu tes the emerging radiation 



199lh for 



211 field for Titan limb rays using the correlated-A; approximation fiLacis and Oinasl . 

212 rapid calculation of atmospheric opacity. The three types of opacity included are: gas vi- 

213 brational bands, collision-induced absorption (CIA) of N2, CII4 and II2, and aerosol/haze 

214 opacity. The ^-coefficients for gaseous bands and CIA opacities are computed and pretab- 

215 ulated over the full range of pressures and temperatures required for the Titan model 

216 atmosphere, as well as the CIRS full spectral range, to enable rapid calculations. The 

217 spectral opacities are pre-convolved with a Hamming apodization function of FWHM 

218 0.48 cm~^, and 50 gf-ordinates are used to approximate the ^-distributions in each spec- 

219 tral bin. 



Spectral line data are taken from the HITRAN 200 4 (for CH4 1304 cm-\ HCN z/2 



221 712 cm"^ and C2H4 z/7 949 cnri 



222 cm , C3H8 V2& band; 



Rothman et a 



Jacquinet-Husson et al 



2OO5I) and GEISA 2003 (for C2H2 z/g 729 



2OO5I ) at lases. We also used supplern ental 



223 sources of line data fo r the C^Hf; z/3 band (822 cm ^ iVander Auwera et al. 



224 i^CHgD (z/6 1156 cm-^ 



Bezard et al. 



20071 ). for 



2003) and for the 1 300-1400 cm"/ ban ds of C3H 



Flaud et al. 



20011 ) that post- 



225 (mostly due to the z/ig 1376 cm~^ and the z/19 1338 cm~^ 

226 date the most recently available GEISA and HITRAN atlases. We made additional scaling 

227 corrections to the line strengths for the 1^26 band of CsHg, and al so to the line width and 



228 temperature dependence parameters for both the GEISA and the 

229 lists, justified in detail in Appendix B. 



Flaud et al. 



mm line 



230 Additionally, for the (1378 cm ^) and z/g (1468 cm ^) bands of C2H6, which are not 
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231 included in any current line atlases, we used a 'pseudo-linelist' made publicly available by 

232 the Jet Propulsion Laboratory (http:// mark4sun.jpl.nasa.gov/pseudo.html). This list is 

233 an approximation which seeks to model the overall absorption of these bands by using a 

234 uniform distribution of 'pseudo' lines of varying strengths, without actual knowledge of 

235 the real line positions or intensities. 



236 To fit the continuum we used two further sources of opacity; (i) co 

237 opaci ty of six molecular pairs of H2, N2 and CH4 (see detailed references in 



lision-inducec 



Teanby et al. 



238 l2006l ). and (ii) a simple haze with a constant (scaleable) number of particles per gram 

239 of atmosphere above the tropopause. The haze had either (A) a uniform grey spectral 

240 absorption, or (B) a spectral absorption cross-section derived by applying Mie theory to 

241 0.2 fim radius spheres (no scattering), using the real and imaginar y refra ctive index data 



242 for laboratory haze analogue (tholin) as measured by 



Khare et al 



fll984f ) 



243 Previously, most m o dels of CIRS 



Vinatier et al 



2007a 



Nixon et al. 



itan mid-infrared limb spectra flTeanbv et al 



2008b 



Jennings et al. 



2008 



Coustenis et al 



2007 



20081) 



245 have used the approximation of a single limb ray originating from the center of the detec- 

246 tor fields of view (FOVs) projected onto the limb: the 'infinitesimal' FOV approximation. 



247 However, 



Teanby and Irwin 



(120071 ) have investigated the spatial smoothing effect of the 

248 small but finite mid-IR detector FOVs when applied to Titan limb sounding. They con- 

249 elude that approximating a limb view using an infinitesimal ray introduces non- negligible 

250 errors when the highest spectral resolution of CIRS (0.5 cm~^) is used and many spectra 

251 are co-added to lower the NESR. 



252 More recently 



Nixon et al 



(I2OO9I ) have analyzed laboratory and in-flight field-of-view 

253 response measurements of CIRS, and computed real 1-D spatial response functions for 

254 all 20 mid-infrared detectors. Using these functions, they investigated systematic errors 

255 introduced by either the infinitesimal ray approximation or the finite (or 'boxcar') ap- 
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256 proximation for Titan limb viewing, relative to using the real detector spatial responses. 

257 Their conclusion was that, for limb averages of as few as 20 spectra, the errors due to 

258 both approximations became significant compared to the noise level. 

259 Therefore, in this work we have included the actual 1-D spatial responses for each 

260 detector as follows. For the n^^ individual spectrum, recorded by detector m, let a„ be 

261 the tangent altitude (in km) at the FOV center projected onto Titan's limb and dn be the 

262 spacecraft-tangent point distance in km. Then the n*^^ spatial weightin g function quantize d 



263 onto a vertical grid Zi of 1-km altitude increments is then (c.f. Eq. 10 of lNixon et al.U2009l ): 



J2iRm{{Zi - an)/dn) 



265 Where Rm is the 1-D (F-aver aged) response func tion of detector m expressed in angular 



266 units (tabulated in Table 4 of iNixon et al.l . 120091 ). For the average of all spectra then 



267 the spatial weighting function is: 

1 ^ 

268 W{Zi) = J-^Yl '^n{Zi) (3) 

n=l 



270 This 'effective' spatial weighting function is used to compute a mean model spectrum, by 

271 weighting synthetic spectra calculated for limb rays at each altitude: 

/(z>)= J2 W{z,)I{z,,u) (4) 

273 where z^nin and ^max are the lowest and highest altitudes where W{zi) ^ 0, and v is 

274 wavenumber. 

275 In practice, the 1 km resolution calculation grid is both unnecessarily fine to achieve 
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276 sufficient modelling accuracy, and also prohibitively slow, as it necessitates the calculation 

277 of ~ 50 individual limb rays to produce a single spectrum at each modelling iteration. We 

278 performed trial retrievals at various grid spacings, finding that 5 km grid point spacing 

279 is sufficient for > 99% accuracy on both FP3 and FP4 (at which level the accuracy 

280 becomes dominated by other effects, such as knowledge of limb pointing), and we therefore 

281 interpolated the function W{zi) at 5 km increments relative to the mean limb tangent 

282 altitude of the average spectrum, cutting off at a minimum weighting of 0.001. 5 km 

283 resolution is approximately 10% of the scale height in Titan's lower stratosphere. 

284 The resampled W{zi) functions in 5 km steps for the FP3 and FP4 spectral averages 

285 analyzed herein are shown in Fig. [2l 

Fig. 2. Appears Here. 



286 3.3 Retrieval Algorithm 



287 The retrieval algorithm is a mathem atical formalism used to find the 'optimal estimate 



288 set of model parameters, by iteration (IRodgers 



20001 ). At each step of the retrieval, the 



289 forward model is used to calculate a synthetic spectrum, as described above, which is 

290 then compared to the data spectrum. A cost function is calculated which is similar to 

291 a test of the difference between the data and model spectrum, with an additional 

292 term to provide smoothing of the solution, also known as the a priori constraint. Model 

293 parameters are then adjusted along the downhill gradients that tend to decrease the cost 

294 function, so that divergence between the model and data will be reduced on the next step. 

295 When a preset convergence limit is reached - i.e. when the percentage change in the cost 

296 function fall s below a threshold value, the retrieval cea s es. Th is algorithm is identical to 

297 that used in Nixon et al. (2008a) - see also Irwin et al. 



(120081 ) for further details. 
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298 4 Results 



Fig. 3. Appears Here. 



299 In this study we first used the 1225-1325 cm~^ spectral region of the z/4 methane band 

300 to retrieve atmospheric temperatures, holding the methane abundance constant. Fig. [3] 

301 shows the result. Note that the retrieved temperatures in the main information region 

302 of the lower stratosphere are significantly colder than the HASI profile, by as much as 

303 10 K - much greater than expected due to latitude variation alone (HASI = 10°S; CIRS 

304 = 2°N). This discrepa ncy has been noted and rema rked on in several previous studies of 



305 CIRS data (see 52 of 



Flasar and Achterberg 



20091 . and references therein), and can be 



306 seen to be part of a larger difference where the HASI stratopause is considerably lower 

307 (by ~ 60km) than that determined by fitting CIRS data. The spectroscopic parameters 

308 for the methane z/4 band are well-studied, and the CIRS data cannot be reconciled with 

309 the HASI temperature profile, so a solution to this problem has not yet been reached. 



The error on the temperature profile includes both spectral noise and also t 



311 unc ertainty on the Huy 



See 



Nixon et al. 



gens GCMS methane abundance quoted by 



Niemann et al. 



ne 5% 



fl2005h . 



313 error variation. 



(j2008al ) for description of the error propagation. Fig. H] shows the vertical 
Fig. 4. Appears Here. 



314 The temperature profile was then fixed, and we carried out numerous separate re- 

315 trievals on specific spectral regions suspected to harbor propane emission bands. In these 

316 retrievals, the grey 'aerosol' opacity was scaled to obtain a fit to the continuum level, and 

317 simultaneously relevant minor gas VMR profiles for each spectral region were adjusted to 

318 fit the emission bands. 

319 Two separate models were calculated for each spectral region: (i) a best fit spectrum 

320 without propane, to reveal the propane signature in the residual, and (ii) for spectral 
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321 regions where CsHg line data were available, we fitted a model for all gases including 

322 propane, to obtain a CsHg abundance, allowing the values from multiple CsHg bands to 

323 be intercompared. The results are now discussed in detail. 

324 Propane Band Detections 

325 Figs. [5] and M show the results of modeling various CIRS spectral regions that were 

326 candidates for propane emissions. In each figure, the left hand column shows the data 

327 spectrum and model (grey haze used), while the right hand column shows the residual 

328 after subtraction of model from data, along with the corresponding portion of the PNNL 

329 laboratory absorption spectrum. The laboratory spectrum has been convolved with a 0.48 

330 cm~^ Hamming kernel for direct comparison to CIRS data, and arbitrary scaling is used. 

Fig. 5. Appears Here. 
Fig. 6. Appears Here. 

331 Four strong bands are clearly detected on FP3 and four on FP4, as indicated on the 

332 figures, along with the locations of several weaker undetected bands. Note that the band at 

333 1053 cm~^ is detected on both FP3 and FP4 where the spectral ranges overlap. The band 

334 assignments for all vibrational modes of CaHg are listed in Table II. At this juncture, we 

335 are obliged to discuss a confusion regarding propane band numbering in the pre-existing 

336 literature. 

337 [TABLE II] 

338 Historically, there has been some inconsistency over the actual band designations for 

339 propane's infrared spectrum. This appears to be traceable to two original studies - a 1965 

340 paper by Gayles and King, hereafter GK65, and a 1972 paper by Shimanouchi, hereafter 

341 S72 - that introduced differing numbering schemes (see Table II). The GK65 scheme was 
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342 apparently followed by lMaguire et al.l (jl98ll ). iGiver et al.l (119841 ) and all subsequent anal- 

343 yses of the IRIS dataset, and thenc e propagated into the Cassini CIRS community. The 



344 S72 scheme however was adopted by 



Flaud et al 



(I2OOII) 



Roe et al. 



(120031 ) and also by the 



345 U.S. National Institute of Standards and Technology (NIST), and appears preferable now 

346 both because of that official endorsement, and also because all 27 predicted vibrational 

347 modes are accounted for. Therefore, we use the S72/NIST numbering scheme in this work. 

348 Table II lists both sets of band designations as an aid to translating band numbers from 

349 earlier works (especially Voyager papers). 

350 Under the NIST numbering scheme, we clearly detect the z/26 (748 cm^^), 1^21 (922 

351 cm^^) and z/20 (1053 cm~^) bands on FP3, which all have strong c entral Q-branch peaks 



198l|). 



352 All three of these bands were claimed to be detected by Voyager (iMaguire et al. 

353 Note that the shape of the 1^21 band appears somewhat different at room temperature 

354 (PNNL spectrum) compared to the colder Titan spectrum, due to the presence of hot 

355 bands on the short-wavelength side of the fundamental in the lab spectrum. The band 

356 on the other hand lacks a Q-branch, and the detection is uncertain at this time, although 

357 we do appear to see a minimum in the right location near the band center at 869 cm~^. 

358 Note that the residual of the ethane fitting appears surpris ingly large, given that we ar e 



20071). 



359 using very recently updated line parameters for the z/g band (iVander Auwera et al. 

360 This appears to be caused by trying to minimize the overall fit? the presence of the 

361 non-fiat continuum. 

362 We must be careful here to emphasize possible aerosol effects: a non-grey aerosol could 

363 easily account for some or all of the gradient in the residual. So long as we cannot positively 

364 distinguish between the two opacity sources, the propane us detection should be regarded 

365 as still somewhat tentative. A clear way to make progress would be to derive spectroscopic 

366 parameters for the band from lab data and include it in the gas model: any remaining 

367 residual gradient would then likely be due to aerosol alone. 
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368 On FP4, the strongest propane feature inside the clean spectral range (1100-1400 

369 cm~^) is the Q-branch of z/ig at 1376 cm~^. The z/7 band (1158 cm~^) is similar in shape 

370 to the Us of FP3, and is also evident here after the modeling/removal of CH3D and ^^CHsD 

371 emissions. In previous works, this band has sor netimes been modeled w ith a band opacity 

372 model to improve the D/H determination (e.g. 



Coustenis et al. 



1989a). Other weaker Q- 



373 branch peaks of z/25 (1192 cm ^) and z/19 (1338 cm ^) fall inside the FP4 nominal range, 

374 but are at the noise level. 

375 Proceeding with caution outside the nominal range (see Appendix A), we see a clear 

376 signal of the Q-branch of z/20 (1053 cm~^), with improved signal-to- noise (S/N) over 

377 the FP3 detection. (We can discount aliasing here, as there is no noticeable emission at 

378 the aliased frequency of 987 cm~^.) In the short-wavelength non-clean band, we also 

379 apparently detect the Q-branch peak of z/24 (1472 cm~^), although it is apparent here that 

380 significant residual remains after modeling which is not attributable to propane. There 

381 may be several factors at work here to undermine our model: (i) the ethane pseudo-linelist 

382 is likely to be imperfect; (ii) the ethane (1468 cm~^) and uu (1469 cm~^) bands remain 

383 strong out to ~ 1600 cm~^, and therefore will contribute significantly through aliasing 

384 ('folding' around 1529 cm~^); (iii) the continuum may have haze opacity that is difficult 

385 to assess in the presence of the other uncertainties. 



386 4-^ Propane Gas Abundance 



387 Having demonstrated the detection of multiple propane bands by CIRS, we wished 

388 to model those for which line data is available. As stated earlier, only the 1^26 band at 

389 748 cm~^ has been previously analyzed to obtain abundances of CsHg on Titan, due to 

390 its unique availability in an existing atlas (GEISA). However, due to more recent work 



by 



Flaud et al. 



(120011 ). line parameters for propane bands in range 1300-1500 cm ^ were 
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392 available to us, allowing the new possibility of retrieving propane abundances from these 

393 bands also. After some experimentation, we concluded that accurate modeling of the 

394 z/24 band at 1472 cm~^ was not possible, due to the uncertainties in ethane and haze 

395 opacity, combined with the complexitity of folding the (unknown) aliased spectrum above 

396 1529 cm~^. However, fitting the z/ig at 1376 cm~^ proved to be a more tractable problem. 

397 We used two spectral windows for the abundance measurements: 740-760 cm~^ {^26) 

398 and 1370-1390 cm~^ ('^is); choosing relatively narrow spectral regions to avoid the effects 

399 of, respectively, the very strong acetylene band at 729 cm~^ and large instrumental noise 

400 spike at 768 cm~^ (FP3), and a non-fiat continuum that showed apparent broad haze 

401 features at the scale of ~20 cm^^ (FP4) . Haze spectral features are the intended subject 

402 of a future study, and are not pursued further herein. 

403 We then added propane gas back into the model atmosphere and separately retrieved 

404 its best-fit abundance for each window using the previously determined temperature pro- 

405 file. Our objective was not to determine the absolute abundance of propane at a single 

406 altitude to high accuracy, where we could have fitted all regions simultaneously, but rather 

407 to compare the relative results for FP3 and FP4. Early retrieval attempts using these win- 

408 dows alerted us to a problem with the GEISA line strengths for the z/26 band, which is 

409 decribed in detail in Appendix B. This problem was remedied by rescaling the spectral 

410 lines prior to the final retrievals. 

Fig. 7. Appears Here. 

411 The final results of the fitting are shown in Fig. [71 The red line shows the residual for 

412 a forward-model calculation using the retrieved abundances of all gases except propane, 

413 which is removed to show its contribution. The other lines show the residual after fitting 

414 with all gases including propane and either grey haze (blue line) or tholin haze (green 

415 line). The fits with the two haze models are indistinguishable over the narrow 20 cm~^ 
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416 window. 



417 [TABLE III] 

418 The individual abundances determined from each model fit are given in Table III, 

419 along with errors. Included in the error calculation (added in quadrature) are (i) the 1- 

420 (J random noise on the spectrum; (ii) propagated temperature retrieval error from two 

421 sources, (a) the spectral noise from the temperature retrieval, and (b) the 5% uncertainty 



in the methane abundance quoted by 



3.5 of , Nixon et al 



Niemann et al 



(120051 ) ■ combined as described in 



(l2008al ): (iii) the uncertainty in the line intensities of propane in the 



424 spectral line lists, estimated at ~10%. 

425 The results show that the z/ig yields systematically higher values than the i^26) com- 

426 patible at the 2-0" error level but not at 1-cr. We suggest that two remaining systematic 

427 effects are the sources of the difference: (i) any residual structure in the continuum at a 

428 scale of less than 20 cm^^ for either spectral region that is not included in our slowly- 

429 varying haze models, and (ii) any imperfections in the pseudo-linelist for the the ethane 

430 z/g band at 1378 cm~^. An actual linelist for the z/g and z/g bands would be very valuable 

431 in eliminating this source of error. 

432 We also note that the z / qr resu. 



433 surement of 



Vinatier et al 



t is in very good agreement with the CIRS limb mea- 
(l2007af ) for 15°S, 3 mbar: ~ 4 ± 1.5 x 10~^, also derived from 



the z/^fi band us i ng a d ifferent (line-by-line) model and retrieval algorithm. Note that while 



Vinatier et al. 



(l2007al ) did not scale the GEISA line intensities by the 0.42 factor used 

436 by us, they partly compensated for the known problem with the z/26 band intensity by 

437 omitting the vibrational part of the partition function when calculating the line strength. 

438 This would tend to underestimate the line strengths, providing a similar correction factor 

439 to ours, and the agreement between the results is not anomalous. 
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440 5 Summary and Conclusions 



441 We have modeled a low-latitude limb spectral average of Titan, removing all gas emis- 

442 sions except propane. By examining the residual, we demonstrate that CIRS has detected 

443 at least six separate bands of CsHg, including the well-known 748 cm^^ band (z/26)- We 

444 also show clearly the structure of the z/7 band at 1158 cm~^ that underlies the band 

445 of CH3D; the z/21 and z/20 bands at 922 and 1053 that were weakly detected by IRIS; 

446 and at least two others in the range 1300-1500 cm~^ (z^ig and z/24) that are evident after 

447 the removal of methane emission. Detection of the weak centered on 869 cm~^ that 

448 underlies the wing of the C2H6 z/g band is uncertain. 



449 We have modeled the 1^26 band 

450 spectral line data from the work of 



using GEISA line data, and also the uiq band using 



Flaud et al. 



(120011 ). None of these determinations are 



451 problem-free: the z/26 band lacks spectroscopic data for the hot band lines at 749.5 cm"^, 

452 whereas the z/ig band suffers from a lack of accurate continuum shape, and possible in- 

453 accuracies due to the use of the JPL pseudo-linelist for the uq band of C2H6, for which 

454 no formal lists are available. Despite these uncertainties, the values agree at a 2-a er- 

455 ror level. We suspect that systematic uncertainties (haze, ethane line list) are still of 

456 greater concern for FP4 than FP3, and therefore we recommend that the lower value 

457 (4.2 ± 0.5) X 10~^, which agrees with previous analysis, is considered valid at the 2 mbar 

458 level for low latitudes probed by these retrievals. 



459 Further work is required to accurately characterize the broad spectral 



"eatures of Ti- 



tan^ mid-infrared continuum, as has been performed for the far-infrared (Ide Kok et al. 



461 120071 ). Additional laboratory spectroscopy of alkanes is also required to support modeling 

462 of planetary spectra. It is vital that high-quality laboratory line lists are prepared for all 

463 bands of propane in the range 700-1300 cm~^, and also the uq and z/g bands of ethane. 

464 The present lack of accurate spectral information is a serious obstacle, not only to accu- 
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465 rate measurement the abundance of propane on Titan and the giant planets, and deriving 

466 underlying haze opacities, but also because our inability to model and remove propane's 

467 prolific signature in many infrared spectral regions is a major hindrance to the search for 

468 the signatures of new stratospheric trace species. 
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A Spectral range of CIRS FP4 



487 The subject of the useable spectral range of CIRS FP4 is of considerable interest in this 

488 paper, due to the existence of many propane bands at the extremes of FP4 sensitivity. 



489 Further details o 

490 report 



Brasunad (120041 ). 



the CIRS spectral bandpass are available in an internal CIRS team 



491 The nominal spectral range is 1100-1400 cm , which significantly underestimates 
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492 the original bandpass: from 1020-1529 cm~^, based on initial Nyquist sampling of the 

493 interferogram signal in the instrument. Frequencies from outside this range are aliased 

494 into the range, according to the formula: 

495 ^'nominal + ^aliased = 2 X i^bandpass i-^-^) 

496 A second-stage numerical filter is then applied that tapers the response towards the 

497 ends of the range, and leads to a definition of the 'clean' bandpass where the strength of 

498 the ahased signal has a relative response strength of less than 1%: 

"99 i^clean = (2 X i^bandpass) ~ ^Ipc (A- 2) 

500 where i>ipc is defined as the (aliased) frequency where the response of the numerical filter 

501 R — Ffp4(^') has dropped to 1% of the peak response: i>ipc = -^FP4~''^(-Rpeak x 0.01). 

502 By this definition, the clean bandpass of FP4 is 1131-1433 cm~^, close to the canonical 

503 wavenumber range (1100-1400 cm~^). 

504 Between the Nyquist and 'clean' bandpass limits at the low wavenumber side (1020- 

505 1131 cm~^), the aliased signal comes from the region 889-1020 cm~^, where there is little 

506 emission except from the relatively weak lines of C2H4. In addition, roll-off of detector re- 

507 sponse at long wavelengths means that the 'clean' bandpass is somewhat underestimated. 

508 These effects mean that the range 1020-1131 cm~^is very useful, especially as the aliased 

509 spectral range 889-1020 cm~^ is known from the FP3 spectrum and its importance can 

510 therefore be directly assessed. 

511 The non-clean range from 1433-1529 is more problematic, because the detector reponse 

512 remains strong, and more critically because we have no knowledge of the signal in the 

513 aliasing range above 1529 cm~^, except from predictations of emission bands of known 

514 gases (e.g. the weak U2 band of methane at 1534 cm~^). Therefore, interpretation of this 
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515 spectral range must proceed with caution. 



516 B Modifications to propane spectral line lists 



517 In this work we used propane line data from two sources: (i) the GE ISA atlas, for the z/26 



Flaud et al.l (1200 ll ) for the bands in the 



518 band at 748 cm^^, and (ii) recent measurements by 

519 range 1300-1500 cm~^. During preliminary intercomparisons of the propane abundance 

520 derived from the two spectral regions, we found large differences in the retrieved values 

521 amounting to a factor ~ 2 (FP4 higher than FP3). We investigated possible causes of 

522 this discrepancy, and discovered that the propane line strengths that are in the GEISA 

523 2003 atlas are incorrectly scaled: this problem is expected to be fixed in the forthcoming 

524 GEISA 2009 release, that should be available at the time of publication of this article. 

525 We also suggest modifications to the Lorentzian Half- Width at Half Maximum (HWHM) 

526 and temperature exponent of the spectral lines. A detailed discussion follows. 



527 B. 1 Radiative Transfer Treatment of Spectral Lines 



528 For atmospheric radiative transfer applications we wish to compute the spectral depen- 

529 dence of the absorption co-efficient fc^ for each molecular energy transition (i.e. spectral 

530 line). Absorption due to all individual lines in a given spectral region is then included 



531 in the overall opac ity calculation, a. 

532 opacity). Following 

533 line is: 



■o ng w ith other sources (particles, collision-induced 



Goody and Yung (119891 ) . the spectral absorption due to an individual 



Sf{u - z/o) 



(B.l) 
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535 where v is frequency, z/q is the unperturbed frequency, / is a normahzed function describing 

536 the hne shape, and S is the hne strength defined as: 



S 



kydv 



:b.2) 



538 The hneshape function may be given by the classical Doppler (thermal) or Lorentzian 

539 (pressure/collision) profiles, or more generally a convolution of both, known as the Voigt 

540 profile: 



{(z/ - z/q) - (nz/o/c)} +01,2 



KI-nkT 



exp 



-mu 
2kT 



du 



(B.3) 



542 In the above formula k, h and c are the usual constants of Boltzmann, Planck and speed 

543 of light; m is the molecular mass, T is temperature and ai is the Lorentzian HWHM. 

544 Therefore, in order to compute k^ we must know the following: m and T, which are 

545 specified by the atmospheric problem; and ai and S, which are parameters obtained from 

546 experimental measurement and tabulated in a database. Both of these parameters have 

547 a temperature dependence, and in addition ai is a function of pressure also. First, the 

548 dependence of az, on atmospheric conditions is specified by: 



To 
T 



:b.4) 



550 where {aL,o, PqiTq) is a laboratory (usually) room temperature measurement, and n is a 

551 free scaling parameter known as the Temperature Dependence of Width (TDW). Second 



552 the line strength as a function of temperature is given by (e.g. 

553 p9033): 



Lacis and Oinaj . 



1991 



^ c: ^0^0 / he 



1 



1 
T 



(B.5) 
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555 where E" is the energy of the lower state of the transition, and R and V are the rotational 

556 and vibrational partition functions (sums over microphysical state probabilities). The 

557 rotational partition function is R{T) = where x = 1 for linear molecules and x = 3/2 

558 for non-linear molecules, such as propane. The vibrational partition function is usually 

559 approximated by a quantum harmonic oscillator so that: 

560 V{T) = n (1 - exp{-huj/kT))~^' (B.6) 

j 

561 where uj is the frequency of the j*'^ vibrational mode, for each of the 27 mode frequencies 

562 given in Table II, and gj is the degeneracy of each mode, equal to 1 for this gas. 



563 B. 2 GEISA Line Strengths for 



564 The spectroscopic line list provided in the current GEISA 2003 edition is unchanged 

565 since first added in 1991, and consists of ~ 9000 lines assigned by S. Daunt (unpublished) 

566 from room temperature spectra. S. Daunt (private communication) has told us that the 

567 integrated band intensity: 



S 



Sidu 



:b.7) 



569 for ~ 3500 lines for the fundamental mode only (excluding hotbands - transitions from 

570 non-ground st ate energy levels) w as scaled to match the band sum for the entire region 



571 as reported 



Giver et al 



Giver et al. 



(119841 ). which included hot bands. The value of S reported by 
(119841 ) is 4.33 x 10~^^ cm molecule"^, and we have separately estimated from 

573 the PNNL propane spectra a value of 4.27 x 10~^^ cm molecule"^, in good agreement. 

574 However, the integrated band intensity of the GEISA lines is slightly lower (3.76 x 10^^^ 

575 cm molecule"^). We have calculated V{2'dQK) = 2.71, which is approximately the ratio 

576 of S for all transitions to ground state transitions only, and therefore, we estimate that 
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577 the S for the i/se fundamental mode should be (1/2.71) x 4.27 x 10"^^ = 1.58 x lO^^^ 

578 cm molecule"^. The strengths of individual lines in the atlas should therefore by scaled 

579 by the ratio of the current GEISA S to the predicted: i.e. by 1.58/3.76 = 0.420. We have 

580 hence rescaled the GEISA lines for the z/26 band by this factor prior to computation of 

581 the fc-tables, which were used in our retrievals. 



582 B.3 Lorentz Broadening Parameters 



583 In the current GEISA 2003 edition, the Lorentz HWHM is given as 0.08 cm ^ 

584 for the z/26 band. However, the FWHM for the pressure broaden i ng of CsHg by N2 has 



585 been separatel y reported as - 119 c m ^ by 



586 0.146 cm by 



Hillman et al. 



Nadler and Jenningsl (119891 ) (at 296 K) and 



(119921 ) (at 175 K). The latter authors also used these two 

587 measurements at different temperatures to infer a value for the n exponent in Eq. IB. 41 of 

588 0.50, which is lower than the 0.75 in GEISA at present. We have therefore adopted values 

589 of a^^Q = 0.12 cm~^ and n = 0.50 for all the bands of propane available to us, including 

590 the FP4 bands. 
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TABLE I 
Titan Observational Data 



Focal 


Spectral 


Latitude 


Mean 


Altitude 


Mean 


No. of 


NESRt 


Plane 


Range (cm~^) 


Range 


Latitude 


Range (km) 


Alt. (km) 


Spectra^ 


(W cm~^ sr~^/cm~-^) 


FP3 


600-1100 


30°S-30°N 


9.5°N 


100-150 


130 


568 


3 X 10-9 


FP4 


1100-1500 


30°S-30°N 


2.1°N 


100-150 


125 


616 


3 X 10-1° 



^" 2 X 10^ space and 2 x 10^ shutter spectra used in calibration. 
^ NESR = Noise Equivalent Spectral Radiance (1-cr). 



TABLE II 
Propane IR Spectral Bands 

Freq. 1972 1965 Mode Laboratory 

(cm~i) No.t No.^ Type Spectra 



216 
268 
369 



14 
27 
9 



cms FPi 

Torsion 
Torsion 
CCC deform 



IR gas-phase inactive 









CIRS FP3 


748 


26 


21 


CH2 rock 


869 


8 


8 


CC stretch 


922 


21 


16 


CHg rock 


940 


13 




CH3 rock 



Husson et al. (1992) 



gas-phase inactive 



1053 



20 



— CIRS FP3 & FP4 

15 CC stretch 









CIRS FP4 




1158 


7 


7 


CH3 rock 




1192 


25 


20 


CH3 rock 




1278 


12 




CH2 twist 


gas-phase inactive 


1338 


19 


14 


CH2 wag 


Flaud et al. (2001) 


1376 


18 


13 


CH3 s-deform 


Flaud et al. (2001) 


1392 


6 


6 


CH3 s-deform 




1451 


11 




CH3 d-deform 


gas-phase inactive 


1462 


5 


5 


CH2 scissor 




1464 


17 


12 


CH3 d-deform 




1472 


24 


19 


CH3 d-deform 


Flaud et al. (2001) 


1476 


4 


4 


CH3 d-deform 


Flaud et al. (2001) 








CIRS ENDS 














2887 


16 


3 


CH3 s-str 




2887 


3 


11 


CH2 s-str 




2962 


2 


2 


CH3 s-str 




2967 


10 




CH3 d-str 


gas-phase inactive 


2968 


23 


18 


CH2 a-str 




2968 


15 


10 


CH3 d-str 




2973 


22 


17 


CH3 d-str 




2977 


1 


1 


CH3 d-str 





''^Shimanouchi (1972) -tGayles and King (1965) 



TABLE III 



Retrieved Propane VMR (ppm) 



Band 


Grey Haze 


Tholin 


748 cm-^ 
1376 cm-i 


(4.2±0.5) xlO-^ 
(5.7±0.8) xlO-'^ 


(4.3±0.5) xlO-^ 
(5.7±0.8) xlO-'^ 



778 Figure Captions 

779 Fig. [1] Initial gas profiles for Titan spectral model. Upper plot: major isotopic species. 

780 Lower plot: minor isotopic species. 

781 Fig. [2] Plots of W{zi) spatial weighting function for limb spectral averages, resam- 

782 pled onto a 10 km grid relative to the mean tangent altitude, and normalized such that 

783 12iW{zi) = 1.0. Dashed line: FP3. Solid line: FP4. The effective weighting function of 

784 the final mean FP3 and FP4 spectra is different, due to two effects: (i) the different real 

785 detector spatial responses, and (ii) the slightly different distributions of center altitudes 

786 of the spectral sets comprising each average. 

787 Fig. [3] A priori (dashed line) and retrieved temperature profile (solid line, with la error 

788 bars) from Titan low-latitude (30°S-30°N) spectral average from limb data (100-150 km 

789 tangent altitude). The grey shaded box shows the approximate region of real temperature 

790 information in the retrieval, as bounded by the min and max values of FWHM pressure 

791 levels of all temperature contribution functions assessed across the entire wavenumber 

792 range. At higher and lower levels the temperature profiles is smoothly joined to the a 

793 priori. A second weaker information region also exists above ~500 km due to methane 

794 hotband emissions. 

795 Fig. m Comparison of various contributions to the error on retrieved temperatures 

796 from fitting the methane z/4 band. Dotted line: initial a priori error estimate. Dashed line: 

797 retrieval error due to spectral noise. Dot-dashed line: error due to methane abundance 

798 profile uncertainty. The combination of these two (spectral and methane profile errors) 

799 lead to the final vertical temperature error: solid line. 

800 Fig. [5] Identification of propane spectral bands from CIRS Focal Plane 3 (600-1100 cm~^) 

801 Left column: CIRS averaged limb spectra (9°N, 130 km), four selected data ranges (black 

802 line) and synthetic spectral fit (red) with all gases except CaHg in model. Right column: 
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803 residual spectrum after model subtraction (red) compared to PNNL laboratory absorption 

804 spectra (blue) and 1-a NESR (grey). Propane bands are identified by text and dashed 

805 lines. Note: incomplete removal of C2ile band emission is evident in the residual (d)at 

806 840-850 cm-i. 

807 Fig. [6] Identification of propane spectral bands from CIRS Focal Plane 4 (1000- 

808 1500 cm~^). Left column: CIRS averaged limb spectra (2°N, 125 km), four selected data 

809 ranges (black line) and synthetic spectral fit (red) with all gases exept CsHg in model. 

810 Right column: residual spectrum after model subtraction (red) compared to PNNL labo- 

811 ratory absorption spectra (blue) and 1-a NESR (grey). Propane bands are identified by 

812 text and dashed lines. Note: spectral baseline in the range 1440-1480 cm^^ (panel (h)) 

813 is not well modeled by grey haze, leading to over-estimation of the ethane abundance 

814 (features marked with '*') as the retrieval attempts to minimize the overall of the fit. 

815 Fig. [7] Residuals of fitting the two propane spectral regions. Red lines: residual of 

816 forward calculation with no propane in model atmosphere, all other gases at previously 

817 retrieved values. Blue line: residual after fitting with propane included, and grey haze 

818 model. Green line: residual after fitting with propane included, and tholin particle haze 

819 model. Grey bands are the 1-a NESR. (a) 1^26 band region, fitted by GEISA 2003 propane 

820 line data modified as described in Appendix B. Note that the propane spectral data does 

821 not include the hot band at 749.5 cm~^, and also the significant difference to the 'baseline' 

822 level due to propane, (b) The i^ig band fitted with the recent spectral data of Flaud et al. 

823 (2001). 
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(a) Data and Model Fit 



(b) CIRS Residual and Lab Spectrurr 
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(a) Data and Model Fit 



(b) CIRS Residual and Lab Spectrurr 
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Nixon et al. Figure [71 
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